Instanton Glass Generated by Noise in a Josephson-Junction Array 
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We compute correlation function of a superconducting order parameter in a continuous model of 
a two-dimensional Josephson-junction array in the presence of a weak Gaussian noise. When the 
Josephson coupling is large compared to the charging energy the correlations in the Euclidian space 
decay exponentially at low temperatures regardless of the strength of the noise. We interpret such 
a state as a collective Cooper-pair insulator and argue that it resembles properties of disordered 
superconducting films. 
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Systems consisting of densely packed metallic grains 
have been studied for decades, see, e.g., for review Ref. 
[l|. They exhibit peculiar electronic properties that stem 
from the quantum tunneling of electrons between the 
grains. Numerous models of Josephson-junction arrays 
have been employed to describe properties of granular 
superconductors 0. The most recent theoretical re- 
search in this area has been inspired by the experimen- 
tal evidence of a sharp low-temperature superconductor- 
insulator transition in disordered films 

Ui, Hi] (see 

also earlier experimental works Q)- Various theoreti- 
cal scenarios of the effect have been proposed: Suppres- 
sion of the superconducting order parameter by disorder 
in two dimensions 8] , Bose condensation of vortices [9| , 
trapping of Cooper pairs due to Coulomb blockade [lO|, 
collective superinsulator phase [llj , insulating Josephson 
phase quenched by the magnetic field [l^ , overheating 
of electrons due to inefficient electron-phonon processes 
[l3l] . and Anderson localization of Cooper pairs 14 1. 

In this Letter we are concerned with the correlations 
of the order parameter in a two-dimensional array of 
strongly coupled superconducting grains in the presence 
of a weak Gaussian noise. We find that regardless of the 
strength of the Josephson coupling between the grains, 
the imaginary-time correlation function of the order pa- 
rameter decays exponentially in 2-1-1 dimensions. We call 
such a state an instanton glass and interpret it as a col- 
lective Cooper-pair insulator in which Cooper pairs are 
weakly localized within areas that include a large number 
of grains. 

We consider superconducting grains that are suffi- 
ciently large so that the fluctuations of the magnitude 
of the complex order parameter A can be ignored. This 
condition is satisfied when the distance between electron 
energy levels in the grain, S, is small compared to |A|. 
However, the phase 9i of A on each grain is a dynamical 
variable described by the Hamiltonian 
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The sum is over all grains, with {i < j) denoting the sum- 
mation over the nearest- neighbor pairs of grains. The 
first term in Eq. ^ corresponds to the charging effect 
due to the Cooper-pair exchange between the grains, with 
Ui being the number operator for the excess Cooper pairs 
at the i-th grain and Eq being the charging energy of the 
grain. The second term describes the Josephson coupling 
of strength Ej between the grains. The form of the 
last term in Eq. ([T]) implies the existence of additional 
weak links that allow some leakage of the Cooper pairs 
in and out of the superconducting grains. In our model 
these weak links are external to the two-dimensional 
Josephson-junction array of the grains. They exist on 
top of the strong links between the grains. We assume 
random distribution of phases and call such a disorder 
"Josephson noise" . 

The prevailing view is that the ground state of a gran- 
ular superconductor depends on the ratio of the Joseph- 
son coupling energy and the charging energy. If this ra- 
tio is large, the Cooper pairs move freely between the 
grains and the system is a superconductor. If the ratio 
is small, which should be expected for small grains with 
large Ec — (2e)^/(2C), then moving an excess Cooper 
pair into the grain costs too much energy. In such a case 
the Cooper pairs are localized on individual grains and 
the system is an insulator. In our model we require 



Ej:$>Ec, 
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so that our granular film is electrically close to a ho- 
mogeneous film. Quantitatively, this condition trans- 
lates into a large dimensionless tunneling conductance, 
g = 2Trh/{e^R) ^ 1, with R being the normal resis- 
tance of the Josephson contact between the grains. In 
this case the charging energy becomes renormalized by 
the Coulomb screening of the excess charge on the grain, 
so that [H 
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Large Ej, provides small differences between Oi and Oj 
at the neighboring grains. On the contrary, we assume 
the Josephson noise to be weak, 



(4) 



so that the corresponding small tunneling conductances, 
allow an arbitrarily large difference between Oi and 
In what follows we treat 4>i as a dynamical random field. 
We show that contrary to the previous findings, in the 
presence of the noise, the Cooper pairs are localized at 
T = within areas of size r cx l/(ej) even under the 
condition 

Since rii and 9i are canonically conjugated variables, 
we have 



ih- 
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This allows one to replace ni in Eq. ([T]) by Mi/(2i?^). 
Then the first term in Eq. JT]) acquires the form of the 
"kinetic energy" . The Euclidean action corresponding to 
the Hamiltonian ([1]) is 
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where r = it. Without the kinetic term this action is 
equivalent to the XY spin model in a random field that 
has been intensively studied in the past [l^. Without 
the last term Eq. ([7|) has been also intensively studied 
(with and without dissipation) in 1980s in connection 
with the possibility of the low-temperature re-entrant su- 
perconductor - normal metal transition due to quantum 
fluctuations of the phase [3]. Superconductor- insulator 
transition at Ec ~ 2£^j in a two-dimensional Josephson- 
junction array has been confirmed by Monte-Carlo simu- 
lations 17| . We are not aware of any theoretical investi- 
gation of the ground state of the model described by Eq. 
d?]) under the condition 

The essential features of the model can be studied by 
considering a square Josephson-junction array with a lat- 
tice spacing a, E"^ — Ec, Ej — Ej at T 0. Small 
difference of the phase for the neighboring grains allows 
one to write for the nearest neighbors 

cos[0(ri) - 0{r,)] « 1 - i[0(ri) - ^(r,)]^ (8) 



(ri) + (r, - r,) • [V0(r)]r=r. 
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Substitution of these equations into Eq. ([T]), summation 
over the four nearest neighbors in the square lattice, and 



replacement of the summation over i by the integration 
according to J2i ~^ I d^f/a^, yields a continuous field 
model described by the action 



dr / d^r ■ 
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It is convenient to use dimensionless variables: 

- X - y - T 
x = - , y = -, T = — 

a a tq 



with 
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In terms of these variables Eq. (jlOp acquires a simple 
form: 



+ ^[l-cos((? -</.)] 
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Here the integration is over dimensionless Euclidian co- 
ordinates {d^r — dxdydf), V6' is the 3d gradient of 9 
with respect to these coordinates. 

We are interested in the limit of Ej ^ Eq when Se/f 
is large compared to h and the phase 9{r) is a well-defined 
semiclassical field. Quantum dynamics of such a field is 
dominated by the extremal trajectories of Eq. (fT3)) satis- 
fying 
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At £,/ = this equation possesses a solution = const 
that describes a global superconducting current. In gen- 
eral, for such a current to exist, the phases at distant 
points must be correlated. We, therefore, want to com- 
pute the correlation function 



C(fl,f2) 

i/gi[e(fi)-e(f2)] 



(A(fi)At(r2)) 



) = (cos[0(fi)-0(r2)]), (15) 



where r = {x,y,f). The average in Eq. is over all 
pairs of points (ri,f2) in 2-1-1 dimensions that are sepa- 
rated by the same distance |fi — f2|. 

The non-linear dynamics of the field expressed by Eq. 
(|14|) usually presents a problem for the computation of 
the correlation function in Eq. (fT5|l . Below we use a 
mathematical trick that under the condition ^ allows 
one to obtain C(f i, f2) exactly with a conventional choice 
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for the noise. We introduce a random two-component 
vector field, 

f(f) = [fiJ-z] - [e,/(r)cos(^(r),ej(r)sin<^(r)], (16) 
and write Eq. (jl4p in the integral form: 
1 



0(f) 
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is the Green function of the 3d Laplace equation, satisfy- 
ing V'^G{f) = S{f). Substituting Eq. ^ into Eq. (USD 
one obtains 

C(fi,f2) = (expj^ J dMG(ri - r) - G(f2 - f )] 
X [/i(r)sin(?(r) - ^(r) cos 0(r)]}) , (19) 

To proceed with the calculation of the space average in 
Eq. (|19p one needs to choose the model of the noise. The 
simplest choice corresponds to the Gaussian distribution 
for the probability, P, of any given realization f (f): 



P[f (r)] cx exp 
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which also provides definition of (ej). With this assump- 
tion Eq. becomes 



C(fi,r2) 
D^{{} exp 
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where / D'^{f{r)} = / D{fi}D{f2} denotes functional 
integration over the realizations of disorder. 

At first glance the evaluation of the above correlator 
may seem hopeless because it requires the explicit knowl- 
edge of 9{r) created by f(f). To see that G(fi,f2) can 
be calculated exactly in the limit of weak noise, we first 
notice that according to Eq. ([4]) and Eq. (fT4|) the contri- 
bution of f (r) to the spatial derivatives of 6'(f) is generally 
small. In fact, as is shown below, the significant change 
in 9{r) occurs over the distances f ^ Ej/{ej) ^ 1. This 
means that the value of at a certain point f has very lit- 
tle correlation with f at that point. Consequently, to the 
lowest order on the noise, the dependence of 9{r) on f (r) 
in Eq. (PT|) can be neglected and the remaining Gaussian 
integration over f can be easily performed. As a result, 



sin0(r) and cos0(r) in the exponent nicely combine into 
sin^ 6{f) + cos^ 9{f) — 1, yielding 

G = exp 1 dV-[G(fi - f) - G(f2 - f)]2| 
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Further integration gives 



G(ri,r2) = exp 



I 



1='-^. (23) 



Phase correlations that decay exponentially in both 
space and imaginary time indicate that the Josephson- 
junction array is in the insulating state [l8| . Since the 
instanton solutions of Eq. describe tunneling trajec- 
tories of Cooper pairs, we call such a state an "instanton 
glass" . The presence of a finite correlation length in the 
(x, y) plane, 



Rc 



(24) 



implies that at any given moment of time the phase cor- 
relation is lost over spatial distances greater than Rc- At 
any point in space there is also a finite correlation length 
in the imaginary time. 
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that implies the existence of the energy gap. 
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characteristic of an insulator [iSj . It represents the local- 
ization energy of Cooper pairs within overlapping areas 
of size Rc- 

Applied to a granular film, the above results mean that 
Cooper pairs are localized within regions of size Rc that 
are large compared to the average size of the grain a. If 
the dimensions of the film, L, are greater than R^, the 
film should be a Cooper-pair insulator. At T ^ A/g 
the conductivity of such a film must be due to the ther- 
mal hopping of Cooper pairs between regions of size Rc, 
obeying the law exp(— A/g/T). So far we have not in- 
cluded the magnetic field into the problem. Its treatment 
within our model is much more involved but a plausible 
speculation can be made about the expected effect of the 
field. The latter is known to suppress Josephson tunnel- 
ing. The weaker the coupling the smaller is the critical 
field that destroys tunneling. In our model the weakest 
links are the external ones characterized by the coupling 
strength ej. They will be destroyed by the field first. 
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making the localization length ([24)1 infinite and destroy- 
ing the insulating gap (j26p . The non-monotonic depen- 
dence of the Josephson coupling on the field [l^l may lead 
to the non-monotonic field dependence of the gap. The 
expected temperature and field behavior of the instanton 
glass is, therefore, in a qualitative agreement with exper- 
imental findings in disordered superconducting films Q. 

Also in agreement with experiment is the fact that such 
a behavior is pertinent to a two-dimensional system and 
does not appear in a three-dimensional superconducting 
sample. Indeed, for a 3d sample, the calculation similar 
to the one performed above, but with the Green function 
G4(f) = -l/(47r2|rn instead of Gair) = -l/(47r|r|), 
provides a power-law decay of the correlations of the 
phase in a 3+1 Euclidian space. The corresponding 3+1 
correlation length is, therefore, infinite and the gap is 
zero. Note that such a dependence of the correlations on 
the dimensionality of space is typical for models with a 
continuous order parameter in the presence of quenched 
disorder [2^,|2H. Our approach is, in effect, an extension 
of the Larkin-Imry-Ma theorem to the Euclidean space- 
time for problems that involve tunneling in the presence 
of noise. 

In conclusion, we have demonstrated that Joseph- 
son noise, no matter how weak, destroys long-range 
space-time correlations of the order parameter in a two- 
dimensional Josephson-junction array regardless of the 
tunneling conductance. We have argued that this ef- 
fect can be relevant to the observed low-temperature 
superconductor-insulator transition in disordered films. 
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